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O ■ Abstract 

Low-dimensional copper oxide lattices naturally manifest electronic states with strong short 
range quantum entanglement, which are known to lead to remarkable emergent material properties. 
However the nanometer scale many-body wavefunction is challenging to measure or manipulate 
in a simple way. In this study. X-ray induced dd electronic transitions are used to suppress spin 
entanglement across a single lattice site in the spin-1/2 antiferromagnetic chain compound SrCu02, 
revealing a class of cuprate magnetic excitations that result from breaking the spin chain. These 
measurements are the first to employ two closely spaced X-ray resonances of copper (M2 and M3) as 
a form of natural 2-slit interferometer to distinguish between different types of electronic transition 
and resolve how they influence the dynamics of nearby spin-entangled electrons. 



Low-dimensional copper oxide lattices naturally manifest strong short range electronic 
interactions that cause remarkable emergent properties, such as high critical temperature 
superconductivity (high Tc) and the appearance of exotic particles defined by the collective 
dynamics of multiple electrons l|-|6|. Directly observing the microscopic quantum entangle- 
ment of these materials requires a specially selected probe that can single out an individ- 
ual atom and perform precise manipulations on a femtosecond time scale to measure the 
many-body collective excitations that occur when the quantum entanglement of neighboring 
electrons is broken or altered. To achieve this goal, we have used resonant inelastic X-ray 
scattering (RIXS) to generate electronic transitions confined to a single copper lattice site 
that disconnect magnetic coupling in the antiferromagnetic spin ^ chains of the quasi-one 
dimensional cuprate SrCu02. Our measurements make use of the closely spaced M2 and M3 
resonances of copper as an analogue to the 2-slit interferometer, and use the pattern of scat- 
tering interference in conjunction with dramatically improved energy resolution to explore 
a previously unobserved class of cuprate magnetic excitations that result from breaking the 
spin chain. Patterns of quantum interference between M2 and M3 resonance channels are 
shown to unambiguously reveal the atomic orbital component and scattering phase difference 
of low energy excitations, critical pieces of information that would otherwise require ultrafast 
coherent measurements of a type not yet developed. Our results demonstrate that similar 
measurements have the potential to reveal quantum interference and magnetic dynamics in 
quantum materials based on Ni, Co, Mn and Fe. 

The cuprate SrCu02 (SCO) has been of interest in recent years as a model realization 
of a one dimensional (ID) spin | antiferromagnet [gHuJ], and as part of a class of quasi- 
ID strongly correlated materials manifesting extremely large nonlinear optical effects 12 |. 
The copper oxide lattice is composed of chains of corner sharing copper oxygen plaquettes 
(Fig. IC), with a | spin moment in the singly occupied 3dx-2-y2 orbital on each Cu atom. 
The chains are formed in pairs that are weakly connected through magnetically frustrated 
edge-sharing plaquette boundaries. Spin interactions along the chain axis are unusually 



strong for a cuprate, with a Heisenberg coupling 



neutron and photoemission spectroscopies [6 



;erm of J/f=0.23eV measured from optical. 



13|. Like many High-T^; superconductors. 



the crystal structure of SCO lacks apical oxygens above and below the Cu site, leading to a 
large crystal field splitting effect differentiating the energies of Cu 3d orbitals. 

The orbital arrangement of electrons surrounding a copper atom can be changed by 



inducing a crystal field excitation. To do this, the energy of incident X-rays is set near the 
copper M2 and M3 resonances, corresponding to the excitation of electrons from spin-orbit 
coupled Cu 3pi/2 and 3p3/2 core levels to the 3da,2_j^2 orbital. Crystal field dd excitations form 
when an electron from a different 3d orbital (e.g. 3d^y in the Fig. IB diagram) transitions 
to the vacated 3p core level, leaving the hole in a new d-orbital as labeled on the RIXS 
spectrum in Fig. lA. The electron added to the 3dx2-y2 orbital is decoupled from the upper 
Hubbard band due to the attraction of the valence hole |l4l . [l5|, causing the excitations 
to be essentially immobile across the Brillouin zone (Fig. lA, inset). Our high resolution 
data reveal that the onset edge of the lowest energy excitation closely resembles a ISOmeV 
Lorentzian (Fig. lA), suggesting that the excited state has a long lifetime of ra.j^~4.3fs 
(/i/150meV=4.3fs). Due to the strong spin interactions of SCO, these dd excitations create 
an unusual situation in which the electronic configuration surrounding a single copper atom 
is reshaped into a metastable state that recovers more slowly than the characteristic intersite 
propagation of magnetic collective modes 6|. 

Since RIXS spectroscopy at transition metal M2 and M3 edges has not been extensively 
utilized, exploratory measurements of the dd excitations across a wide range of incident 
photon energies are presented in Fig. 2A (bottom). Summing the RIXS intensity within 
1-3. 5eV excitation energy results in a curve representing the total probability that a dd mode 
of any symmetry will be excited, as a function of incident energy (Fig. 2 A, top). Two humps 
centered on the M2 and M3 resonance energies can be fitted using closely spaced Lorentzian 
functions with a 1:2 intensity ratio as expected from the relative core level degeneracies. In 
contrast, separately tracing the resonant scattering intensity of individual dd modes in Fig. 
2E reveals distinctively different patterns. The Sd-s^^-r^ dd excitation primarily resonates in 
a single peak around hz/=76eV, while the 3dxy excitation resonates over a broader energy 
range with two distinct peaks approximately 3eV apart, more widely spaced than the actual 
energy gap between 3pi/2 and 3p3/2 core levels. 

The considerable overlap between M2 and M3 resonance suggests that the dd modes are 
excited via intermediate states in which the core hole occupies a quantum superposition of 
the 3p i /2 a nd 3p3/2 levels. Within the widely successful mean field "direct RIXS" scattering 
model 15l4l9l| , overlapping the resonant scattering signal from closely spaced core levels will 
produce an effect similar to conventional two-slit interference. The interference pattern is 
not caused by wavefronts intersecting in space, but rather by the phase of quantum paths 



through 3pi/2 and 3p3/2 states in the energy/time domain (Fig. 2D) as presented in the 
following scattering equation: 



R{hu, dd) oc 






hu - Em2 + iT/2 hu - Em3 + «r/2 

Resonant scattering intensity of a selected excitation mode {Kihu ,dd)) is deter- 
mined based on different real-valued quantum amplitudes and phases for scatter- 
ing from M2 and M3 (e.g. ^M2,xy, 0M2,xy for the Sdxy mode), as well as 
the resonance energies (EM2=74.8eV,EM3=77.55eV) and inverse resonance lifetime 
(r=2.5±0.1eV=/i/(0.26±0.3fs)). A wide range of information about the M-edge scatter- 
ing process can be determined from these fitting parameters which will be discussed further 
in the Supplementary Online Material (SOM), including the Ms-resonance per-atom scat- 
tering cross section which we find to be larger than for the higher energy L3 resonance by 
a factor of roughly two. As shown in Fig. 2E, the quantum interference effect described 
by Equation (1) nicely reproduces the squeezing together and spreading apart of resonant 
intensity observed in the data. 

The distinctive forms of these resonance patterns occur because the spin and orbital 
transitions involved in creating a dd excitation through the M2 and M3 intermediate states 
differ by a phase that shifts the constant phase contours of light scattered from M2 and 
M3 along the time axis (arrows in Fig. 2D). Using the resonance energetics parameters 
itted from the data, one can predict the full resonance profile in the direct RIXS process 
15| and obtain scattering phases of m2-G mz^t^ for '^^■iz'^-r'^ excitations and ~0 for Sd^..^^ in 
agreement with those fitted in Fig. 2E. Adopting a lowest order approximation that the 
dd excitations are Lorentzians along the energy loss axis, the resulting numerically derived 
spectrum is a good qualitative match for the experimental data, as can be seen in Fig. 
2A where the two are overlaid. Previous studies have noted the existence of these classes 
of quantum interference, but have not established or made use of them experimentally by 
tracing the incident energy dependence of individual excitation modes [171, l20[. In these 
measurements we experimentally observe that the phase shifts picked up during resonant 
scattering provide a new element of spectroscopic information to determine the orbital and 
spin symmetry of RIXS features which can otherwise be a subject of controversy and model- 



dependent speculation |l5l . 
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In RIXS experiments performed with resonance at more tightly bound core levels (e.g. K- 
edge), higher order "indirect RIXS" processes involving a large continuum of highly complex 
intermediate resonance states are known to create patterns of incident energy dependence 
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that are not found in M- resonance data 16|, ll7| • In indirect RIXS processes, 
the dynamics excited by a scattered X-ray are biased by the choice of incident photon energy 
and can become difficult to relate to desired correlation functions. These mechanisms cannot 
be responsible for the resonance patterns observed in our M- resonance data (Fig. 2E), as 
the characteristic indirect RIXS charge transfer excitations of SCO are not present (see 



B 



Fig. 4C and Ref. (8|). The phase interferometric resonance patterns identified in our 
data are much simpler, caused by just two intermediate resonance states that have the same 
configuration of low energy electrons regardless of the incident energy used for measurement. 
Because of this, the dynamics induced by each separate dd excitation are almost entirely 
independent of the incident photon energy, and the quantum interference pattern can be 
used as a tool to selectively enhance a specific excitation mode that one wishes to measure 
without compromising the excitation line shape. 

Of the excitations seen in our data, the "id^y dd mode has the longest lifetime and is 
expected to cause the greatest disturbance to the many body electronic structure of sur- 
rounding lattice sites. Calculating the resonant phase interference pattern reveals that the 
Sdjjy excitation can be observed with maximal intensity in the YZ scattering plane and is 
uniquely enhanced relative to other features at X-ray energies of hi^~74.5eV, just below the 
M3 resonance (Fig. 3A). At this photon energy, we rotate the sample azimuthal angle by 90" 
to selectively suppress the Sd^j; excitation relative to other features, so that its line shape 
is observed 'in isolation' by taking the difference between measurements in the YZ and XY 
scattering planes (Fig. 3C). The spectra measured in these two geometries can be directly 
compared because M-resonance RIXS involves very limited momentum transfer from low 
energy photons (Q~0) and a nearly featureless absorption cross-section dominated by a flat 
non- resonant background (see SOM). The isolated "idxy dd mode has an asymmetric line 
shape, also somewhat visible in our lower resolution K-pre-edge data (Fig. 3C, inset), with 
a sharp rising edge and shallow trailing edge that cannot be attributed to the charge and 
spin features known from earlier RIXS studies 15l. Il7|. 

The line shape of dd excitations encodes the rich many body dynamics brought about by 
disturbing the complex electronic structure of a cuprate. However, limited resolution has 



previously made it impossible to distinguish dd mode profiles from the Lorentzian contours 
predicted in single atom mean field models [l5|,|l7|. To understand where the non-Lorentzian 
asymmetry in our 3d., data comes from, we consider the inherent qnantnm flnctuations of 
a spin i Heisenberg antiferromagnet, which are particularly strong in ID materials [ij] such 
as SCO and are proposed to provide the glue for cuprate superconductivity in quasi-2D 
lattices. Introducing a localized "idxy dd excitation to an antiferromagnetic chain in SCO 
breaks spin interactions by doubly-occupying the Cu orbitals that interact with neighbor- 
ing spins (3da;2_y2 and 3d3^2_j,2), causing the final state to be composed of two broken spin 
chains that terminate to the right and left of the excitation site. An exact diagonalization 
simulation of spin dynamics after breaking the SCO spin chain is shown in Fig. 4A, using 



the spin interaction term of J//=0.23eV p, |9|, O, ll3|. At the instant when the excitation is 



created, magnetic moments close to the excitation site are correlated in a short range antifer- 
romagnetic pattern that does not change significantly over the core hole lifetime (5t<0.3fs). 
As time evolves towards the full lifetime of the ?)(\xy dd excitation (r~4.3fs), adjacent spins 
lose most of their correlation with the original spin of the transition site and the magnetic 
configuration up to 6 lattice sites away (2.4nm) changes significantly. Excited magnetic 
harmonics of the simulated broken chain are therefore seen to have significant dynamics 
over the lifetime of a dd excitation, making them important in determining the spectrum 
observed by RIXS. Calculating the 3dxy excitation spectrum within an ideal J//=0.23eV 
spin chain reveals that it closely matches the experimental line shape due to a high energy 
tail of spin excitations (Fig. 4B). 

To evaluate the importance of spin excitations in other features as well, the entire spec- 
trum is fitted using identical parameters with spin dynamics turned on and off in Fig. 
4C(right). When spin dynamics are discarded by setting J//=0, or equivalently by treating 
the system as a classical Ising antiferromagnet, qualitative deviations are observed for all 
three dd features. This clearly shows that excited magnetic collective modes of the broken 
spin chain are expected to contribute significantly to the dd excitation spectrum. These spin 
modes occur due to quantum spin fluctuations and are not a feature of mean field spin mod- 



els, un 
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ie the local spin flip excitations that can effect the spectrum of the 3d322_r2 mode 



181 ]. Similar line shape effects have been observed previously related to attosecond 



timescale charge relaxation near a perturbed site in simple molecules 29|], but magnetic 



dynamics in a cuprate occur on a comparatively slow femtosecond timescale and would be 
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challenging to observe without high energy resolution achieved at the M2- and M3- resonance 
channels. 

Another fundamental effect governing the line shape is the lifetime of RIXS charge exci- 
tations, which has been implicated as causing charge excitation profiles to be consistently 
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23|. Strikingly, our 



broader than 5E>150meV in late transition metal oxides 
high resolution data reveal that creating RIXS excitations on isolated spectator ions such as 
Nd in nickelate Nd2-xSr2;Ni04+5 can lead to a much sharper line shape (AE=110±10meV 
in Fig. 4B, inset). Contrasting these two cases suggests that the intrinsic width observed by 
RIXS may be increased by electronic interactions within the transition metal oxide lattice. 
To consider dynamical modes that dd excitations may interact with, the Ms-resonance dd 
excitation spectrum in Fig. 4C is superimposed with measurements of excitations across the 
copper-oxygen charge transfer gap determined by K-edge RIXS and optical ellipsometry. 
The leading edge of the 3dxy dd excitation does not overlap with a large density of charge 
transfer excitations, and has a sharper line shape than any other feature in the spectrum. As 
the energies of the dd excitations increase beyond the charge transfer gap energy, a broad- 
ening trend is observed. The 3d::i,2 ^2 excitation, which is consistently the broadest dd mode 



across many classes of cuprates 15|, has the greatest quantum overlap with charge transfer 
excitations and can easily decay into a delocalized charge transfer collective mode. 

The overall picture that emerges is one in which dd excitations suppress spin interactions 
on a single lattice site for long enough to change the local correlation and singlet-like entan- 
glement between spins |3[, and then decay due to dynamical interactions with the surround- 
ing crystal. The spin excitations triggered by a dd mode occur solely due to non-classical 
quantum spin fluctuations and are thus expected to become weaker in more classical spin 
systems such as 2D cuprates and materials with larger local magnetic moment [ij]. Some 
of the intensity in the modeled RIXS line shape can be attributed to freely propagating 
double- and quadruple- (etc.) spinon excitations, however other changes in the spin state 
are localized near the dd excitation site and do not necessarily resemble free spinous. 

These spin modes induced by dd crystal field excitations are a new kind of many body 
excitation that can be explored with high energy resolution and excitation tunability from 
quantum interference at the M- resonances, and demonstrate how X-rays can interact with 
a complex quantum system by creating an immobile dd excitation 'plug' that blocks interac 
tions through a single lattice site. Our study also reveals a general method for observ- 



ing the phase component of RIXS excitations by using the closely spaced transition metal 
M2 and M3 channels as a natural interferometer, a measurement that would otherwise re- 
quire coherent phase-resolving experimental techniques that have not yet been implemented. 
Looking ahead, our data and analyses demonstrate a number of distinct features and con- 
crete advantages to M-resonance RIXS over the much better known L- and K- resonance 
RIXS spectroscopies that will be important for future measurements, including an enhanced 
per-atom scattering cross section, strong quantum interference, simpler scattering matrix 
elements and a cleaner spectrum with fewer high-order excitation features from indirect 
RIXS. 

Methods summary: 

High resolution M-edge measurements were performed at the BL4.0.3 (MERLIN) RIXS 
endstation (MERIXS) at the Advanced Light Source (ALS), Lawrence Berkeley National 
Laboratory. The data were recorded by a VLS x-ray emission spectrograph equipped with a 
commercial Andor CCD detector cooled to -75°C. Large single crystal samples were cleaved 
in situ and measured at a pressure of 3x10"^'^ Torr at room temperature. All data were 
obtained within 36 hours of cleavage to minimize possible aging of the cleaved surface. The 
resolution-limited peak width at half maximum of the elastic line is better than 5E^35±2 
meV for all M-edge data sets. High resolution measurements, with energy resolution of ~90 
meV, at the copper K-pre-edge were performed using MERIX spectrometer at the Advanced 
Photon Source IXS 30-ID beamline, and more comprehensive measurements including map- 
ping of the K-edge dipole excitation spectrum were performed at SPring-8 BLll-XU. Nu- 
merical simulations and matrix elements relevant to RIXS scattering are discussed in the 
online SOM. 
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FIG. 1: Massive single-atom excitations: (A) Excitations measured using copper M3 resonant 
scattering on SrCu02 are primarily crystal field transitions (highly localized electron-hole pairs), 
labeled by the 3d orbital symmetry of the hole. Statistical error, similar to the data point size, is 
roughly constant across the spectrum due to dark noise of the detector. An inset shows momentum 
dispersion of the Sd^z/yz excitations measured at the K- pre-edge. (B) Ground, intermediate and 
final state {\g),\i),\f)) electron configurations are shown for the direct RIXS scattering process 
by which dd excitations are created using a 3p3/2 (M3) core level resonance. The energy of the 



3d 



x^-~y^ 



orbital in the absence of a core or valence hole is indicated by a dashed line, showing 



that dd excitations are decoupled from the upper Hubbard band. Intra-orbital interactions of the 
3dri.2_y2 and 3dxy electrons are neglected for representational simplicity. (C) The crystal structure 
of copper-oxygen chains in SrCu02 is shown, with axis labels following the convention of two 
dimensional cuprates. 
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FIG. 2: M2-M3 quantum interference: (A) (bottom, blue) Inelastic excitation spectra are shown 
with incident photons spanning the copper M2 and M3 resonance energies, and overlaid with (red) 
a theoretical prediction for the dd modes labeled, (top) The summed inelastic intensity from each 
RIXS curve is compared with a two-Lorentzian fit. (B) The orbital transitions involved in creating 
a 3dxy excitation through M2 and M3 core hole resonances are shown. (C) The constant phase 
contours in a conventional 2-slit interference experiment. (D) The constant phase contours from 
scattering through M2 and M3 resonance are offset in time by a phase that depends on the final 
excitation symmetry, generating distinctive interference patterns for (left) 3d322_^,2 and (right) 3dxy 
dd modes. (E) (circles) Incident energy dependence of the 3dxy and 3d3^2_r2 modes is fitted from 
Equation (1) (solid line) with phase shifts 6m3,z2 = Gm2,z2 +'''" and 6M3,xy = (^M2,xy A dashed line 
shows the fit result with quantum interference disregarded, and arrows highlight the change from 
quantum interference. 
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FIG. 3: Quantum interference and the excitation line shape: (A) A contour plot of the 
SrCu02 M-resonance scattering data is overlaid with colors indicating the amount by which in- 
tensity is increased or decreased by quantum interference based on the orbital symmetries. An 
optimal energy for studying the 3dxy mode (hz/=74.5eV) is indicated with a black arrow. (B) 
The scattering geometry is shown. (C) Changing the scattering plane by rotating the sample 90° 
around the azimuthal (y-) axis suppresses the intensity of the lowest energy dd transition without 
greatly changing momentum because Qi~Q2~0. An inset contrasts (green, high resolution) the 
shape of a single dd excitation revealed by M3 geometrical contrast with (black) the shape of a dd 
excitation observed at the K- pre-edge (resolution (5E~100meV). 
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FIG. 4: Excitation dynamics in a broken spin chain: (A) Time evolution of magnetic 
correlation with the original spin of the dd excitation site is indicated by dots, and the percentage 
change relative to the t=0 magnetic structure is indicated with color shading. (B) The expected line 
shape of a single dd excitation in the SrCu02 spin chain is obtained by diagonalizing a 20-site ring 
and (blue curve) applying 5E=150nieV lifetime broadening, (green circles) The experimental line 
shape of the 3dxy dd excitation obtained by contrasting the RIXS signal in two different scattering 
geometries is shown with statistical error bars. An inset overlays the SrCu02 excitations with 
(silver) the // RIXS spectrum from Nd spectator atoms in Ndi.67Sro.33Ni04+5. (C) (left) The 
dd transitions of SrCu02 measured with M3 photons (hi/=75eV) are fitted as excitations in a 
spin chain, using the inverse lifetimes (SFi) indicated. Dipole induced charge transfer excitations 
measured from (open circles) K-edge RIXS and (black line) optical ellipsometry from Ref. [l3|] are 
overlaid, revealing delocalized charge excitations that dd crystal field modes can decay into, (right) 
Setting the spin interaction (J) to zero gives a prediction that neglects quantum spin fluctuations. 
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